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Abstract. As contrast to the Maxwellian distribution, equilibrium distribution of 
the second kind or turbulent equilibrium is shown to exist under tertiary molecular 
chaos hypothesis to replace the classical binary chaos by Boltzmann. It is expressed 
as bimodal Maxwellians each mode differing by plus/minus root-mean-squares of 
fluctuations in macroscopic variables. Chemical reaction rates calculated using the 
turbulent-equilibrium are obtained in a closed form, accounting for discrepancy 
between experiments and classical theory based on Arrhenius' law that underestimates 
the burning rate considerably. The key issue is the correct estimation of the high-energy 
tail of the distribution function that is minor in population, yet playing a major role 
for reactions with high activation energy. It is extremely sensitive to turbulence level 
in the temperature, causing slightly subcritical molecules to clear the potential barrier 
to cause reactions just like quantum particles undergoing tunnelling effect owing to the 
uncertainty principle. Variance of the fluctuating turbulent chemical reaction rate is 
also calculated, verifying that relative variance based on the turbulent equilibrium is 
low, whereas its classical counterpart (Arrhenius) is pathologically high. A closed set 
of equations governing reactive turbulent gases is presented on this sound basis. 



1. Introduction — A perspective on the classical and turbulence- 
corrected Boltzmann formalism 

It is well-known that each of the transport processes in fluids, namely, diffusion, heat 
transfer, momentum transfer due to shearing motion, and chemical reaction is enhanced 
by the presense of turbulence. The physical origin that causes such appreciable hike 
in those rates is the fractal nature of turbulence|jl[]: it is ascribed to drastic increase 
in net contact area of adjacent bulk of fluid, taking place through the bounding 
surface penetrating into each other area, thereby making it easier to transport mass. 
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momentum and energy. It results in drag rise of a flat plate several times the laminar 
counterpart, also in flame propagation velocity elevated by a few ten times0 under 
internal combustion engine environment(Fig.l). Obviously the latter is a favorable facet 
of turbulence characteristies, whereas the former is detrimental from the viewpoint 
of transport vehicle technologies. In fact, the evidence that the rate increase due to 
turbulence is much higher for mechanism generating thrust rather than drag tempts us 
to imagine a midget planet where airplanes would have to install disproportionally huge 
engines, under the environment of no turbulence. Thus chemical reaction rates, among 
other transport processes, must be influenced by turbulence most sensitively, so is worth 
looking into most detail. 

In general there are two aspects of approach to this problem: The one is to pursue 
the instantaneous value of a physical quantity W_ which is stochastic and fractal, using 
direct numerical simulation (DNS) for example. In principle, any of these quantities are 
described in terms of microscopic density advocated by KlimontovichQ, which we will 
call hereafter the K-formalism. The DNS currently being used in numerical simulation 
of fluid flows is the macroscopic version of the K-formalism 0. It has been applied to 
problems of combustion science, enabling to compute wiggling premixed flame front 
with sufficient resolvability||]. The state of the art, however, is yet to go far to be able 
to provide information such as turbulent flame speed for engine designers. 

Another approach to turbulent combustion is based on an axiom of statistical 
mechanics that claims equivalence of solving for instantaneous value IK of a stochastic 
quantity and for a set of averaged quantities 

{W, WW', W'W'W', ■ ■ ■} (1) 

In the above W = W_\s the average value, 

W' = 'K-W (2) 

is the instantaneous fluctuation subject to W = 0, and fluctuation correlations WW 
etc., refer to those between different points in the independent variable z, e.g., W = 
W(z), W = W(z), etc. The second formalism is the basis on which to found statistical 
mechanics for example; in fact, if one identifies W_ with the microscopic density, then 
its average W represents Boltzmann's function. So the latter may well be called B- 
formalism as contrast to the K-formalism as defined above. 

A flow of turbulent gas is dually stochastic; the one is microscopic (molecular) level 
with the average taken with respect to the molecular velocity, whereas the other is 
macroscopic (fluid-dynamic) level where the average is taken over a volumelet whose 
size is greater than the Kolmogorov scale yet smaller than fluid-dynamic characteristic 
length. 

If one identifies W with chemical reaction rate, there are three levels of descriptions: 
Following the K-formalism W_ is expressed as a sum of on-off step functions depending 



3 



on whether each molecular collision clears the potential barrier to effect inelastic 
collisions 0. Next level is the classical B- formalism where the molecular average alone is 
employed. Then W is shown to be expressed as Arrhenius' law|^, |^. It is coupled with 
fractal analysis to estimate turbulent flame surface area to evaluate turbulent burning 
rate|T^. The third stage, namely, the B-formalism to be discussed here is intended to 
encompass both molecular and turbulent averaging procedures in its framework. 

A series of earlier works along this line[|TT|, [T^ has shed some light on the 
conflicts between the classical Arrhenius kinetics and experimental data attributable 
to temperature turbulence using the exact solution that represents simplified reality. 
In Fig. 2 is shown shock-tube data of low-temperature (spotty) ignition of hydrogen 
oxygen premixture as compared with the classical theory {6T = 0), also with turbulence- 
corrected one, a primitive form|[ll[] of the present theory. Also is shown in Fig. 3 flame 
velocity of turbulent premixed gas||l4| as compared with existing experiment and with 
renormalization group theory. 

Although the earlier theory has thus enabled to predict some features of turbulent 
combustion as to their dependence on turbulent intensity without relying on any 
empirical parameters, no information has been provided regarding their dependence 
on the scale of turbulence. Revised version in the basic formalism has been proposed 
in ref. [0, with applications to nonreactive turbulence, where full informations on the 
scale of turbulence have been made available. 

The objective of this paper is three-fold: The one is to show the exact solution 
for turbulent equilibrium still valid even without the condition of only temperature 
turbulence prevalent. The other is to remodel the previous theory so as to include scale 
of turbulence into formalism for reactive flows now at issue. The third is to show that the 
B-formalism, an average-oriented formalism, provides a sound description of turbulent 
combustion if based on turbulent equilibrium, which otherwise has been considered as 
out-of-date. 



2. Turbulent equilibrium 

In paralled with the classical kinetic theory providing molecular basis of nonturbulent 
fluid dynamics, molecular description of gases based on the microscopic density 

N 

l{z)=Y.S{z-z,{t)) (3) 

s=l 

has given the ground for the direct numerical simulation of turbulent gases p, 0. In 
the above expression, z = {x, v) is a point in the phase space, namely, physical (x) plus 
molecular velocity (v) space, and S denotes the (6D) delta function. This function gives, 
by deflnition, the instantaneous number density in the phase space, so is identifled with 
unaveraged Boltzmann function. It was flrst applied to the plasma kinetic theory^ 
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where the equation governing / is Vlasov's equation. For monatomic (neutral) ideal 
gases, the governing equation, derived as an equation of continuity in the phase space 
with precise spatial and temporal resolvability of Hamiltonian mechanics level, proves 
to be an unaveraged Boltzmann equation! 



where J is the classical collision integral operator 



f{f'f'-ff)VdQdv\_^ (5) 

J . X=T. 



In the above / = f{z) is the microscopic density of the collision partner molecule z, 
f = f{z') likewise with z' specifying the pre-collision state of the collision partner 
leading to z after the collision, V = \v — v\, and dQ is the differential collision-cross 
section of the molecular encounter. 

The microscopic density, being averaged over a space with a number of unidentifiable 
molecules included, defines the Boltzmann function /; 

/ = Z- (6) 

This is where the statistical concept is introduced for a quantity of deterministic 
mechanics. Similarly, two-point correlation function ilj{z,z) is defined using the same 
averaging concept as 

^P{z,z)=ll-ff (7) 

Higher order correlations can be defined likewise. 

The classical kinetic theory is founded on Eq.(4) averaged over the same space as 
(6) together with the Boltzmann's (binary) molecular chaos hypothesis 

^ = (8) 

that is the Boltzmann equation in the classical sense; 

Under the equilibrium condition {d/dt = 0, d/dx = 0) the equation reduces to 

J{z;z)[ff] = (10) 



whose solution is known as the Maxwellian distribution 

^ = ^°=(2^^"P("^) ^ 
W = V — u 

c = {KTlMf^ 



where n, u, T are the number density, fluid velocity and temperature, TZ and Ai denote 
the universal gas constant and molar weight, respectively. 

If a gas is turbulent, fluid variables exhibit (long-range) correlations, so the classical 
binary chaos (8) ceases to hold. Then Eq.(9) is to be replaced with 

W)=[^t + '^-^)f- ^(^; ^)[// + ^)] = (12) 

Thus the equation is not closed at this level: The equation of next step to govern ip 
is 
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/'5(/) + 5(/)f = (/' = /-/) (13) 

which is shown to be identical with the two-particle hierarchy equation of the BBGKY 
theory [ll7|. The actual form of (13) reads 

(I + + A) ^) = ^^''^ + (14) 

+'ilJuiiz,z,z)] + J{z;z)[f'ilj{z,z) + fipiz, z) + il)iii{z,z,z)] 

The simplest possible truncation of the chain of equations at this level is to put 

All = (15) 

namely, to invoke tertiary molecular chaos to replace the classical binary chaos (8). 

Eq.(14) subject to condition (15) turns out to be separable into those for respective 
variables in terms of the wave number 

^{z, z) = R.F.P J (j){z, k)(j)*{z, k)dk (16) 

where symbol * denotes complex conjugate, / denotes a length characteristic of the flow 
geometry and R.P. stands for the real part. Upon its substitution into (14) we have the 
equation governing as 

iuj{k)(f) = i^o(0) 

id d\ ^ \ (17) 

n,i<p)^i- + v. — U-j{z-,zm + <Pf] 

where Lv{k) is the separation constant having dimension of the frequency, and connected 
with the wave number through phase velocity Vp as 

uj = Vp-k (18) 

iFTom Eqs.(17) and (18) the following relationship is seen to hold; 

0*(fc) = 0(-fc) (19) 
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A remark should be mentioned on a renovation that has brought the scale of 
turbulence into formalism: In the older theory variable separation (16) is effected in 
the domain of frequency uj in place of wave number k. Information on the scale of 
turbulence is missed in this form, which may be retrieved in the following way |l5|] : Here 
turbulence is characterized as a wave with frequency uj and wave number k obeying 
dispersion relationship (18). We may assume a plane wave 

= e*^-^$(0, k) (20) 

and discuss its amplitude If further we introduce Fourier transform 

$(2, k) = — — / e-''^-^F{z, s)ds (21) 



f27r)3 



oo 



to work with 'eddy' space s, and employ separation rule written in terms of wave number 
k, namely after (16), a simple formula follows: 

1 f°° 

M^, z) = — — / F{z, s)F{z, s + r)ds (22) 

where r = x — x is the distance between the two space points. Any informations on scale 
of turbulence can be obtained from this correlation formula by taking fluid moment of 
(22), for example, the velocity- velocity fluctuation correlation as 

u'j{x)u'i{x + r) = {ppY^ j WjWiipudv dv (23) 

where w has been defined by (11). 

It is a well-posed problem to seek the equilibrium solution of the new set of equations 
(12) and (17). It stands for a 'turbulent equilibrium' solution to substitute the local 
Maxwellian for nonturbulent gases. This is to solve simultaneous integral equations 

J{z-z)[ff+ I M*dk] = Q (24) 



J(z;z)[/0 + 0/]=O (25) 

The exact solution under the assumption that turbulence prevails only in the 
temperature has been obtained in a closed form by Sagara|12[ by summing an infinite 
series in the temperature fluctuation 5T = (T'^)^/^. The same result has been reached 
using much simpler method|I^]. 

We can show that the exact solution also exists for realistic turbulence where the 
fluctuations in number density 5n = {n'^Y^'^ and fluid velocity Suj = {u^Y^'^ (no 
summation convention here) are also prevalent. Put 

f = f = \{U + fo) (26) 
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(27) 



where /o is the local MaxweUian defined by (11), and is given as follows; 





n 


(27rc±2)3/2 




= n ± (5n 








= Vj-uf 




— Uj ± 




= R(T±5T) 


i? 


= n/M 



exp 



w 



±2 



2c 



±2 



(28) 



It is easily confirmed that the form of (26) is designed to be consistent with definitions 
of the average quantity (n, w, T), and the form of (27) with that of the variance 
(5n, 5T) of each quantity provided that the following relationship holds; 



f 



KK*dk = 1 



(29) 



Direct substitution of (26) and (27) into (24) and (25) reads, respectively, 
\j{z;z)[Uf^ + /o"/^] = \{J+ + J-) 



where has been defined as 

J^^Jiz;z)[fofo] 



(30) 



(31) 



This integral vanishes in view of the fact that J{z;z)[fofo] = 0, also that the bilinear 
Maxwellians in (31) consist of the same family in fiuid parameters (Q.E.D.). 

The deduction above claims that the state of turbulent equilibrium is represented 
by the Boltzmann function of the form (26), along with correlation function 



fo) 



V'(^,^) = 4(/o+-/o-)(/^ 
or the two-point distribution function 

fu{z,z) = ff + ij{z,z) 

~ Tiifofo + /o /o ) 



(32) 



(33) 



At a glance difference between the two Boltzmann functions fo of (3) and /° of (26) 
looks trivial, because the difference is of the first order smallness in 6. The point at issue. 
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however, is that they differ appreciably in the population of high energy molecules that 
play a crucial role in chemical reactions, particularly when the activation energy is high. 
In Fig. 4 is shown the population density of molecules having absolute molecular velocity 
V with turbulence only in the temperature. The Maxwellian distribution [ST = 0) is 
seen to have lean high energy tail, namely, underestimates population of high energy 
molecules compared to reality of turbulence. 



3. Turbulent chemical reactions 

Chemical reaction rate of elementary reaction of the form 

[a] + [b] — > [c] + [d] ([ ] ; chemical symbol) 

usually obeys Arrhenius' law; 

W = An^n^ exp{-E/nT) (34) 

where ria and denote partial number densities of reactants [a] and [b], A is the 
frequency factor and E is the activation energy per mol. Rate law (34) is warranted 
on the classical kinetic theory using a simple collision model. A most clearcut model 
among various ones proposed in late 40's through 50's is due to Present @ : Let the 
intermolecular potential of a molecule [a] relative to a molecule [b] be repulsive, of a 
volcano shape having a crater with radius vq at altitude e. Let also the locus of the 
[b] molecule relative to the [a] molecule written in axisymmetric polar coordinate be 
[r(t),x(t)] where x is the angle between the symmetry axis parallel to V^a = — v^, 
and the position vector r = — x^, then we have 



(mt/2)(r2 ^ ^2^2) ^ g(^) ^ (mt/2)l^2 



(35) 



relationships representing conservation of angular momentum and of energy, where 
X = dx/dt etc., = (m~^ + m^^)~^ is the reduced mass, b is the impact parameter, 
and e(r) denotes the intermolecular potential. Since the radius of the closest approach 
r = is given by the condition dr/dx = 0, the critical collision that separates elastic 
from inelastic collisions is such that the closest approach occurs on the ridge line of the 
crater ; = Tq. This condition determines the critical impact parameter bc{V) from 
(35) as 
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Instantaneous reaction rates for respective species are expressed in the following 
form, 

f (37) 

J'mel 

Utilizing the following formula 

2e 



(38) 



/inel " "V mW^ , 

and transforming integral variables from (va, ^b) to (C7, V) defined by 

U = (maVa + "^b^^b)/("^a + mb) 
V = Vb - t^a 

we have for (33) 

^ dVv[l-^A l{z,)l{z^)dU (39) 

where 

Vc^{2e/m)y/^ (40) 

denotes the critical relative velocity for a head-on collision. 

Chemical reaction rate in classical sense is obtained from (39), upon averaging and 
employing classical chaos (8), as 



w 



. dVvil-^] J_^f{z,)f{z,,)dU (41) 



(42) 



If the Maxwellian equihbrium (11) is employed Eq.(41) is integrated out to give 

Wo{T, na, nb) = An^n^e-'/'^ 
A = {SkT/nm^f^Tirl 

that is nothing but Arrhenius law (34) with activation energy E and the frequency 
factor A written in microscopic parameters. 

It should be noted that turbulence counterpart of the above deduction faces with 
difficulty if conducted phenomenologically, whereas it is simply straightforward if one 
follows the molecular approach as employed here. In fact, if one assumes that the 
instantaneous reaction rate is Arrhenius form (34) 

W_ = Anji^exp{-E/nT) 

and its average given by 



W = Arijiy, ex.p(-E/nT) (43) 
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Decomposing the r.h.s. term under average into means and fluctuations we have 



where (5 = E/IZT is the Zel'dovich number and r = TZT' / E stands for the temperature 
fluctuation. Leading term in the Taylor expansion of the expression under average sign 
gives for /5 1, 



c?r" n\ ' 

which imphes that convergence of the series is extremely slow since the Zel'dovich 
number for combustion reaction is of O(IO^). For example we need to take 23 terms 
before diminishes to unity for (3 = 10. It also means that we need to know 22 

consecutive self correlations in the temperature fluctuation. Needless to say, therefore, 
that the first term (Arrhenius) or a few-term approximation is far from satisfactory. 

On the other hand, if one adopts microscopic approach, namely, starts with the 
averaged version of (39) insted of (43); 




W = 7rr'J dVV l--f, fJAU (44) 



and employs turbulent equilibrium (33) for f^f^ = /n, a simple calculation leads 
immediately to the final expression 

W' = \{W^ + W^) (45) 

with 

= Wo{T ± 6T, ± 6n,, ± 5ny,) (46) 

where Wq has been defined by (42). Formula (45) for the turbulent reaction rate turns 
out to be a bimodal Arrhenius: It allows us to draw a intuitive picture that a half of the 
chemical reaction takes place under temperature elevated and population of reactant 
molecules denser by their root-mean-square turbulent intensity, and another half under 
those quantities lowered by the same amount. 

The bimodal Arrhenius low (45) is the full turbulence version derived on the basis 



of exact solution of Eqs. (24) and (25), supplementing the previous one|TT|, |T^, [T3|] , 
where turbulence prevails only in the temperature. As already sketched it is obvious 
that anomalous increase in the reaction rate is overwhelmingly due to the temperature 
turbulence that feeds extra high energy molecules capable of reaction (see Fig.4). The 
first term of (45) looks as if the potential barrier is lowered by the factor of (1 + 5T/T)~^; 
a macroscopic equivalent of tunnelling effect of quantum mechanics where the potential 
barrier is lowered by a factor of Planck's constant due to the uncertainty principle. 

On the other hand, fluctuation in the reactant number density does not have such 
effects on the reaction rate beyond the flrst order smallness. However, when coupled 
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with the temperature effect, it leaves with a possible account for below-critical lean 
combustion that actually occurs in the operation of the direct injection engines (GDI). 

Rate law (45) also tells us that turbulence in the fluid velocity has no influence on 
the reaction rate (as it should not). In fact, the collision integral for a pair of molecules 
about to collide are not influenced by the fluctuation in the fluid velocity, because they 
are aboard on the 'same boat' in the turbulent ocean. Its effects on turbulent flame 
velocity are exclusively via other transport processes, namely, turbulent diffusion and 
turbulent heat transfer. (See Table 1, next section). 

4. Reaction rate variance 

It is an experimental evidence that instantaneous chemical reaction rate W_ = W+W at 
a fixed point undergoes fluctuation over a wide range, therefore, quality of a statistical 
theory employing any average concepts hinges crucially on reliable estimate for its 
variance around the average W. 

This issue again is beyond the reach of phenomenologies and we have to address 
to statistical approach as developed in the preceding section. To be remarked at this 
point is that what appear in the governing equations are not the variance itself, but 
the fluctuation correlations Z'W where Z stands for fluid variables such as velocity, 
temperature, etc.. 

The equation governing W is obtained from (39) by separating it into mean and 
fluctuating parts. We have, then. 



Since, in general, fluid variable fluctuation is expressed as a moment of fluctuation of 
the Boltzmann function in the velocity space as 



W = 7i„el [/'(^a)/(^b) + /(^a)/'(^b) + /'(^a)/'(^b) - /'(^a)/'(^b)] (47) 



where we have defined 





(48) 



we have, from (43) and (44), 




(49) 



Let the variable separation be invoked after (16) as 




(50) 
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and 



/a/a/b = R.F.I'' I UK)Mka)Mkb)5{K + ka + k^)dk,dkadk^ 



= R.P.l^ J (/);(fe)(/)a(fc)(/)b(fe - k)dk dk 



where, in the last row, —ka is replaced with fc, also fca with fe, and condition (19) has 
been made use of. Then Eq.(49) leads to 



dk(l)*^{k)\ gw{k) - /inel[0a(fe)/b + /a0b(fc) 

/'CO ^ ^ ^ 

dk<p^{k - k)(pb{k)dk] \ = 
-oo 



(51) 



+R. 



If one decomposes gw and into a plane wave after (20) 

„ _ „ik-Xry 

gw — e (jrvK 

and deals with its amplitude Gw and or its Fourier-transform version 
Gw{x,k) = j—-j_^dse-''^-^qw{x,s) 



(52) 



(27r0 
1 



oo 

oo 



we have two alternatives for the reaction rate fluctuation from (51), 



(53) 



Gw{x, k) = /inel 



$a/b + /a^b + R.P.i' / dk^,{k - k)^i,{k)dk 



qw{x, S) = Iinel[FJb + /a^b + ^a^b] 



(54) 



(55) 



depending on which of the wave number (fe) or the eddy (s) spaces is more convenient 
to work with. 

The actual form of qw associated with turbulent equilibrium (26) and (27) is 
straightforward. In fact, in view of (52) and (53), we have 



with 



F^{z, s) = S{x, sy- 



S{x, s) = (27r0'/' 



oo \ 1/2 

,2 ^ 



go ds 



(56) 
(57) 



where go is related to variance of the density fluctuation by 

1 f°° 



(58) 
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Utilizing (56) in (55) we have 

S 



qwix, k) = - 



where 



= A I (na + 5na) (nb - Sub) exp 
+ (ria - (5na)(rib + Sri},) exp 



-E/n (t + 



-E/n T - 



ma - mb 



5T 



6T 



(59) 



(60) 



ma + mb 

Once we have obtained exphcit expression for qyy variance of the reaction rate can 
be calculated from (see Eq.(50)) 

1 



gwgw dk 



(61) 



where (52), (53) have been taken into account. In Fig. 5a is shown the relative variance 
{W"^y/'^ /Wq of reaction rate around classical Arrhenius, which amounts to quantities 
exceeding 0(1) by far, indicating extraordinary fluctuation intensity to be observed by a 
fixed-point measurement. Also is shown in Fig. 5b the same ratio for turbulent reaction 
rate (41). In this case the relative variance stays within sound realm, and approaches 
unity with increase in turbulence in the temperature and with Zel'dovich number. This 
is a sign of stochastic signals characteristic of on-off type. 



5. Phenomenology-coupled Boltzmann formalism on turbulent 
equilibrium 

Statistical description of stochastic phenomena where relative variance far exceeds unity 
is pathological. It is for this reason that a formalism using the averaged Arrhenius 
law fails in describing turbulent combustion. However, as asserted in the preceding 
section, such formalism using the bimodal law as the averaged reaction rate deserves 
reconsideration. For this purpose it is advisable to follow the recipe prepared for 



nonreactive turbulence ||15||: We start with phenomenological equations generalized to 
reactive gases written in instantaneous quantities: 



Ao 



m^ 



df>_ drrij. 
dt dxr 



PU.r 







(62) 



(63) 



4- 



dm.: d 



+ 



dt dx. 



mjTrij, 



(64) 



14 



p. = -yU 



d fm.^.\^ d fElj\ 2 ^ d / ^ 



9Xj \ P J dXr \ P J 3 "''^SXfc \ p 



(65) 



2 

E = pe + =^, H = E+p 
- 2p 



(66) 



(67) 



d / 



(6^ 



Ma,r = -PU 



- dXr 



(69) 



In the above, p, tti^, p and e denote density, massflux density, pressure and specific 
internal energy (thermal part only), p. , Q and M^r ^^re viscous stress, heat flux density 



-jr' 



and diffusion mass flux density for species a, and p^, and are partial density, 
mass fraction and zero-point enthalpy of species a, respectively. Also, p, A and Z2a are 
coefficients of viscosity, thermal conductivity and diffusion for species a, respectively. 
Note that all the equations have the form in accordance with a rule (rule I) that they 
are written in terms of quantities proportional to the density. This is to use mass-flux 
density instead of fluid velocity, pressure instead of temperature, and internal energy 
per unit of volume instead of speciflc internal energy. 

A few words must be mentioned about modifled Pick's law (69) that is not exact 
consequence of the kinetic theory of gases as distinct from other molecular transport 
relationships. The 'effective' diffusion coefficient Da is expressed in terms of binary 



diffusion coefficient for each pair of constituents as 

a, = (i-n)/5:z};Jrb 

b 

= 1/E^abX' 



(70) 



where 



Fb/(1 — ^a) is mass fraction of species b of the mixture with species a 
excepted. The approximate formula (69) with (70) warrants validity for practical use 
except a minor flaw that condition ^Ma.r = is not exactly met. 
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Equations for the average are provided by 

Aa = I^ = (a ; 0, j, 4, a, b, • • •) (71) 
whose actual forms are given as follows: 

Pjr = (Pj>)ns + pu'jK (74) 

Qr- = (Q^Fourier + P<^' (76) 

Ma,. = (M,,,)Fick + P<>? (78) 

W^a=^(W^oL+W^O-a) (79) 



where pjr, Qr, M^^^ and Wg, are average transports of momentum, thermal energy, 
species mass due to diffusion and chemical reactions, viewed from the center of gravity 
frame of reference, respectively. Of the four transport processes three but chemical 
reactions have a structure such that molecular and turbulent transports are additive, 
whereas these are tightly coupled and built in a single exponential form for chemical 
reactions (Table 1). A few remarks are in order in deriving Eqs.(72) throTigh (78) 
from (71): They are exact to 0{Z"^) if Mach number is small enough. Single term 
turbulence corrections as derived here would not be available unless one followed rule 
I plus the following one (rule II) that fluctuation of each density-proportional quantity 
be written in density-independent quantities except fluctuation in density itself, for 
example, = rrij + pu'j + p'uj + p'u'j — p'u'y It is remarkable that turbulence 
correction to each of Navier-Stokes, Fourier and Pick laws is represented by such a 
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simple expression |]I9| without employing the so-called mass average concept. Note also 
that artifice of this kind is not necessary to reach the same result if phenomenologies 
are discarded and the whole deduction to fiuid equations starts from the B-formalism 
in the phase space, which, unfortunately, is practicable only to nonreactive monatomic 
gases. 

Equations governing turbulent transports to close the system emerge from the 
following equations [P^ I 



Z'^K'p + K'^Z'p = Q, ((«, /?); (0,j, 4, a, b, ■■■)) (80) 

where A'^ = A^, — Aq,, and Z'^ denotes fiuctuation of fiuid variable Za. Regardless of 
the fact that those equations are nonlinear, they are again separable exactly in the 
form of (22) and (23), resulting in the following set of equations (see Appendix for its 
derivation); 

Dqo + drqr = (81) 

Dqj + dr{ulqj + q^o^rj + qrj) + qrOu^/dXr - p~^{dp/dxj)qo = 
qjr = -n[dj{p-^qr) + drip-^qj) - {2/3)Sjrdk{p~^qk)] 

-[dul/dxj + du]/dxr - i2/3)Sjrdul/dxk]-^^ 

+p-%qr - {Sjr/3)ql] 



— ^Dqio + dr ( —^q^oul + qrkk ) + H e° q^^w = 
7-1 V7-1 /a 

qrkk = -i\/R)[dr{p~^q4) +p''^q4^dRT/dXr] 

+ [7/(7-l)](i?Tg, + p-Vg4) 

Dq^ + dr{ulq^ + ga,r) + qrdY^/dXr - qwa + P'^qoW^ = 

+ p'^qs,qr 



(l^,r = -pD^dr[p ga) - 



s 

1W = ^ 



dT pR dxr 
±; species [a] produced/disappeared) 
S 



(83) 



^4) 



g4o = g4 + pRqo 



i5) 



where 

dj = d/dxj + d/dsj 
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Dqa = dqa/dt - Vprdq^/dsr 

Of these equations terms of fluctuations in transport processes are listed in Table 2. 

Turbulent transports having appeared in equations for the average are written in 
terms of g's as 

1 



pu'jU'j. 



ph'u'^ 



p(27r)3 
1 



qjqy ds 



7 



p(27r)3 7 - 1 7_oo 
1 



q^^qr ds 



p(27r) 
1 



q^qr ds 



^ ' pR 



^4 dsj 



X 1/2 



(86) 



Thus the two sets of equations for the average (Eqs.(71)) and separated equations for 
fluctuation-correlations (Eqs.(80)) are coupled through those quantities. They have to 
be solved simultaneously subject to homogeneous boundary conditions for g's as 

(a 7^ 40) 



qa[Xh, S) 

{dq^o/dxn)x=x^ = 
at solid boundaries and at laminar-turbulent flow boundaries x — x-^,, also 

qa{x, ±oo) = 

a necessary condition to secure convergence of the integral in the s-space. 



(87) 



(88) 



6. Concluding remarks 

So far the 'Reynolds average' concept is considered not to be applicable to turbulent 
combustion gasdynamics. This negative view has its origin in poor predictivity for 
one of the four average turbulent transports, namely, the chemical reaction rate. It is 
concluded that if one employs bimodal Arrhenius law to replace the classical one, sound 
'average' description of turbulent combustion may well be expected despite its on-off 
structure of signals to be observed at the reaction front, a sign of marginal variance 
amplitude allowed for statistical description. 

The bimodal law as the average turbulent reaction rate and its variance are derived 
on the basis of the turbulent equilibrium, an exact solution of non-equilibrium statistical 
mechanics. It is only through the microscopic approach that the correct estimate for 
population of high energy (reactive) molecules is possible, and thereby high reactivity 
under turbulent environment is elucidated. 

Also proposed are Reynolds averaged gasdynamic equations for turbulent 
combustion with those renovations incorporated. They are coupled with another set 
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governing turbulent fluctuations through the turbulent transport processes to constitute 
a closed system. 
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Appendix. Equations governing turbulent fluctuation-correlations 

The actual expression for of Eq.(80) subject to rules I and II are 
A'o = dp'/dt + {d/dxr){p4 + (Ju\ + p'< - = 



= p = {pUr + p'u'^)l P 



(Al) 



d d 



+p'UjUl + pu'jU'j. — pUjU'j. +p'5jr + (Pjr)Ns\ = 



(Pjr)NS 



—p[du'^/dxj + du'j/dxr — {26jr/S)du'f^/dxk] 
-{dp/dT)T'[dul/dxj + du]/dxr - {26jr/3)dul/dxk] ^ 



(A2) 



. / 1 dp' d 

= 7^ + ^ 

7 — 1 OT ox 



1 



{p'ul + pu'^ + p'u'^ — p'u'^) 



7-1 

+ (Q;)Fourierl+EeK = 
-' a 

(<5;)Fourier = -X&T' / dXr - {dX/ dT)T'dT / dXr 



(A3) 



Al 



dpY: , d 



+ 



pulYl + - pF>; + (M^,)Fick 



dt dx. 



Fick 



^^dYj d{pD.) ^,dY, 



' dXr 



dT 



dXr 



(A4) 



Eqs.(80) are again separable into respective independent variables x and x with full 
nonlinear terms retained when subjected to the following separation rules: 



dk gc,{k)g*p{k) 

/OO POO ^ ^ ^ 

dkg}{k) / dkg^{k-k)g^{k) 
-OO J — OO 

Z'^Z'f^Z'^ = / dkg^{k) / dkglik - k)g;{k) 

J —OO J —OO 



(A5) 



This separation rule preserves the property that periodic part in gf's be separated out 
after (52), namely, 

g^{x, k) = e^^-^G^{x, k) 
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which converts Eqs.(78) in 6D space (a;, x) into those in another 6D space {x, k) as 

D{k)Go + dr{k)Gr = (A6) 



D{k)Gj + dr{k){ulGj + G^oSrj + Grj) 

+(du]/dxr)Gr - p-\dp/dxj)Go = 
Grj = -p[dj{k){p-'Gr)+dr{k){p-^Gj) - (2/3)5,>9,(fc) (p-^Gfe)] 



— [dul/dxj + duj/dxr — {2/3)6 jrdul/dxk] 
+{i/p)[r{GjGr) - {Sjr/S)r{G,Gk)] 



G4 dp, 
'pRdT 



(A7) 



^D{k)G,o + dr{k) i^^G,oul + Grkk \ +J2elG,w 



7 



.7 



dX 



Grkk = -(X/R)[dr{k)(p-'G4) +p-'G,—dRT/dxr] 
+[7/(7 - l)][RTGr + p-'r{GrG,)] 



(A8) 



Co, 



D{k)G, + dr{k){ulG, + + GrdY,/dxr - Gw. + — Wa = 

p 

G.,r = -pDMk)(p-'G^) - ^^^^ + p-'nC^Gr) 

dl pK OXr 

where we have defined the following quantities, 
D{k)G = dG/dt-ik-VpG 

dr{k) = d/dxr + ikr 

T{G^G^) = / Ga{k - k')Gp{k')dk' 



(A9) 



(AlO) 
(All) 

(A12) 



Gao — Ga ~\~ RTGn 



(A13) 



Eqs.(A4) through (A7) govern the 'wave' functions Gq, Gj, (^4(^40) and Ga_, Gy,, 



from which turbulent correlations Z'^Z'f^ are calculated. Relationships between the two 
groups of variables are the following; 

/ Go \ 



p 

pu'j 
pRT' 

P' 
pYi 



Gn — 



G40 

\ Gwa J 



(A14) 
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Turbulent transports that appear in the equations for the average are 
pu'jU'^ = -R.P.l^ / GjG*dk, (Reynolds stress) 

P J —oo 



1 7 



P7- 1 

A- a r p 



/oo 
GiPldk, (turbulent heat flux density) 
-oo 



GaG*dk, (turbulent diffusion) 

X 1/2 



1 / roo \ 1/2 

(r2)V2= /3/ GiGldk) , (turbulent chemical reaction) 

pR \ J -oo J 



>(A15) 



pR 

These equations as derived above in the fe-space are integro-differential equations 
having nonlinear integrals F. These integral expressions are of convolution type, so are 
eliminated by Fourier transform 



Ga{x, k) 



dse~'^'^qa{x, s) 



(A16) 



{2nl)^ J-oo 

leading to an equivalent set of equations in eddy (s) space simply through a 
transformation rule posted in Table A. Their actual forms are: 



Dqq + drQr = 

Dqj + dr{ulqj + g4o'^rj + qrj) + qrduj/dxr - p^^{dp/dxj)qo = 
Qjr = -(J^[djip~^qr) + dr{p~^qj) - {2/2,)5jrdk{p~^qk)] 
-[dul/dxj + du]/dxr - {'2/^)5jrdu\/dxk]^^ 



(A17) 



(A18) 



7-1 



-0^40 + dr 



qAQUl + qrkk 



.7-1 



d\ 



Qrkk = -{X/R)[dr{p ^q^) +P ^qA-^dRT/dXr 
+ [7/(7-l)](i?Tg, + p-Vg4) 



(A19) 



-D^a + driulqe, + ^a.r) + QrOY^/dXr - qWa. + P ^QoWa. = 



-pD^dr{p ^ga) - 



d{pD^) q4^ (9Fa , „_i 



S 



qm = 54 + pRqo 



dT pR dxr 
±; species [a] produced /disappeared) 

S 



(A20) 



(A21) 
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Figure captions 



Fig.l Dependence of flame velocity ut/ul on turbulent intensity u' /ul of methane- air 
premixture in elevated pressure environments. Kobayashi et al (1996) [^. 

Fig. 2 Failure of the classical theory {5T = 0) in predicting ignition time r at 
low-temperatures (T < 1000°i^'), as corrected by turbulent reaction rate formula 

{6T > o)|n . 



Fig. 3 Turbulent flame velocity ut/ul as dependent on turbulent intensity u'/ul of 9.0 
percent if2-air premixture due to a predecessor theory of the current one as compared 



with existing experiments, also with Arrhenius law equivalent and renormalization group 
theory. 

Fig. 4 Maxwellian equilibrium /q {6T = 0) and Turbulent equilibrium /° {6T > 0) 
distribution function plotted against absolute molecular velocity, revealing appreciable 
difference in the population of high-energy molecules eligible for chemical reactions. 

Fig. 5 Relative variance of reaction rates under turbulence around (a) averaged 
Arrhenius reaction rate and (b) turbulent (bimodal) reaction rate. 



Table captions 



Table 1 Molecular and turbulent transport processes 
Table 2 Fluctuations in transport processes 

Table A Conversion rules from wave-number space (fc) to eddy space (s) 
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momentum (shear) 
flux density 


(Navier-Stokes' law) (Reynolds stress) 
[u] = rrij/p) 


heatflux density 


qj- A^^^ +pu'^h' 

(Fourier's law) (turbulent heatflux density) 
(Tt = p/pR) 


partial mass flux 
density 


(Pick's law) (turbulent diffusion) 
= PJP) 


chemical reaction 
rate 


(turb 


A(na + 5n^) (rib + 5nb) exp ^ ^ j 

M(na 5na)(nb 5nb) exp ^ 7e(r-(5r))] 
ulence-corrected Arrhenius law) 



Table 1 
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fluctuations in 


molecular 


turbulent 


momentum (shear) 
flux density 


-{2/3)Sjrdk{p-'qk)] 
— [dul / dxj + duj 1 dxr 


+p-\qjqr - ((^j>/3)gfe] 


heatflux fiensitv 


-{\/ R)[dr{p-\,) 
dX 

+p-'q4^dRT/dxr] 


+ [7/(7-l)](i?Tg,+p-ig,g4) 


partial mass flux 
density 


-pD^dr{p~^qs.) 
dT pR dxr 


+p~^q^qr 


chemical reaction 
rate 


f - + ^{w^ + - wr)] 



Table 2 



wave number (k) space 


-> eddy (s) space 


Ga{x, k) 


qa{x, s) 


ikj 


d/dsj 




-> (9„- — cilcixA + 3 /3s A 


D{k)Ga = dGa/dt - ik-VpGa - 


Dqa^ dqa/dt -Vprdqa/dsr 


/ Ga{k')Gp{k - k')dk' 





Table A 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0001060vl 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0001060vl 



This figure "fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0001060vl 



This figure "fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0001060vl 



This figure "fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0001060vl 



